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The anisotropy parameter (^2), the second harmonic of the azimuthal particles distribution, has 
been measured with the PHENIX detector in Au+Au collisions at ^/s NN — 200 GeV for identified 
and inclusive charged particle production at central rapidities (\r}\ < 0.35) with respect to the 
reaction plane defined at high rapidities (\r)\ = 3-4). We observe that the vi of mesons falls below 
that of (anti)baryons for pt > 2 GeV/c, in marked contrast to the predictions of a hydrodynamical 
model. A quark coalescence model is also investigated. 

PACS numbers: 25.75.Dw 



Event anisotropy is expected to be sensitive to the 
early stage of ultra-relativistic nuclear collisions at the 
Relativistic Heavy Ion Collider (RHIC). The possible for- 
mation of a quark-gluon plasma (QGP) could affect how 
the initial anisotropy in coordinate space is transferred 



into momentum space in the final state. The anisotropy 
parameter V2 for a selection of produced particles is de- 
rived from the azimuthal distribution of those particles. 

dN 

— oc 1 + 2v 2 cos2(</>-$ RP ) (1) 
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where <f> is the azimuthal direction of the particle and 
<I>rp is the direction of the nuclear impact parameter 
("reaction plane") in a given collision. Measurements 
of the parameter v 2 in RHIC collisions have been per- 
formed [1-6] for charged particles and for identified par- 
ticles. The current work reports results for charged par- 
ticles versus transverse momentum (px) out to 5 GeV/c, 
and extends previous measurements for identified parti- 
cles out to 3 GeV/c for n and K, and to 4 GcV/c for 
protons. (Previous measurements of the v 2 for n, K, and 
p extended to 1 GeV/c at ^/s^7 = 130 GcV. [2]) Detailed 
measurements of the azimuthal anisotropy are important 
to eventually discriminate among different possible sce- 
narios for its physical origin. Such scenarios include: hy- 
drodynamical flow of compressed hadronic matter, the 
production of multiple mini-jets, and an anisotropy de- 
veloped during an early quark-matter phase of the col- 
lision. It has been observed that v 2 saturates at px ~ 

2 GeV/c and above [4,5]. The cause of this saturation 
is not yet known; however, we note that at this momen- 
tum the particle composition is very different than at low 
momentum in that the proton yield is comparable to the 
pion yield [7]. This makes the measurement of v 2 for 
separately identified particles especially interesting. 

The measurements described here were carried out in 
the PHENIX experiment at RHIC [8]. About 28 M min- 
imum bias Au+Au collisions at y / s^ = 200 GcV from 
the 2001-2002 run period (Run-2) are used in the anal- 
ysis. Charged particles are measured in the central arm 
spectrometers (|?7| < 0.35) [9] where PHENIX has ex- 
cellent particle identification capabilities [10]. The drift 
chamber (DC) and the first pad chamber plane (PCI) 
together with the collision vertex define the charged par- 
ticle tracks. In order to reduce background, the recon- 
structed tracks are confirmed by requiring matching hits 
in the outer detectors, i.e. the third pad chamber plane 
(PC3) and the electromagnetic calorimeter (EMCAL) or 
the timc-of- flight detector (TOF). In this analysis, the 
TOF detector is used to identify charged particles up to 
4 GcV/c in px- Particle time-of-flight is measured us- 
ing the TOF with respect to the collision time defined 
by beam counters (BBC), and is used to calculate mass 
squared using the particle momentum and the flight path 
length [7]. The timing resolution of the system is ~ 120 
ps. A momentum dependent ±2 a cut on mass squared 
allows particle identification in the px range 0.2 < px < 

3 GeV/c for pions, 0.3 < px < 3 GeV/c for kaons, and 
0.5 < pt < 4 GeV/c for protons. The contamination of 
mis-identified particles is less than 10%. In addition to 
collision time, the BBC provide z-vertex position infor- 
mation. The two beam counters are located at |z|=1.5 m 
from the collision point, covering |?7| = 3 <~ 4. They 
consist of 64 photo-multiplier tubes (PMT) equipped 
with quartz Chercnkov radiators in front surrounding the 
beam pipe. The large charged multiplicity (a few hun- 
dred) in \i]\ = 3 <~ 4 and the non-zero signal of event 



anisotropy in this rj range enables us to estimate the az- 
imuthal angle of the reaction plane in each event using 
the BBC with full azimuthal angle coverage. 

Since the v 2 parameter is in effect a quadrupole mo- 
ment, the anisotropy which gives rise to a non-zero v 2 is 
often referred to as an "elliptic flow." It is extracted by 
first determining the reaction plane angle <I>rp for each 
event, 



tan 2<I>rp 



Sn ch sin 20 PM x 
Sn ch cos 20pmx 
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where n c h is the number of charged particles per PMT 
(determined from the pulse height in each PMT) and 
4>pmt is the azimuthal angle of each PMT. Then, it is 
calculated by the Fourier moment v 2 = (cos 2(<j> — 3>rp)) 
over all particles, for all events in a given sample [11]. 
Corrections [11-14] are applied to account for finite reso- 
lution in the reaction plane determination, and for possi- 
ble azimuthal asymmetries in the reaction plane detector 
response. The bottom-left panel in Fig. 1 shows the av- 
erage cosine of the difference between the two reaction 
planes defined by the beam counters at r\ = 3 <~ 4 and at 
rj = —4 ~ —3 using the the elliptic (second) moment defi- 
nition. In order to improve the reaction plane resolution, 
a combined reaction plane is defined by averaging the re- 
action plane angles obtained from each BBC, using the 
elliptic moment in each. The estimated resolution of the 
combined reaction plane [11], (cos 2($ mcasurod ~ $ truo )), 
has an average of 0.3 over centrality with a maximum of 
about 0.4. The estimated correction factor, which is the 
inverse of the resolution for the combined reaction plane, 
is shown in the top- left panel in Fig. 1. 

The present technique is distinguished by defining the 
reaction plane angle using particles at high rapidity when 
measuring v 2 for particles at mid-rapidity. Other mea- 
surements of v 2 for mid-rapidity particles at RHIC have 
used reaction planes defined with mid-rapidity particles; 
or have employed a technique of measuring angular corre- 
lations between pairs of particles at mid-rapidity. While 
these different approaches generally seek to measure the 
same thing, they are not identical and a variety of physics 
effects can cause them to yield different results from the 
same collision sample [15,16,4]. Because of the large 
rapidity gap between the reaction plane and the mid- 
rapidity acceptance of about 3 units, it is expected that 
this analysis is less affected by non-flow contributions. 
However, we do not observe any substantial difference 
between the v 2 results shown here and published results 
for the v 2 of charged particles at RHIC in the px range 
where they are available. 

The centrality of each collision is defined using the si- 
multaneous measurement of the total number of particles 
measured in the BBC and the total energy measured in 
the zero degree calorimeter [17]. The middle panel in 
Fig. 1 shows the centrality dependence of v 2 for charged 
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particles measured at mid-rapidity (|r/| < 0.35) with re- 
spect to the reaction plane defined above. The centrality 
is measured in percentile from the most central collision. 
The V2 parameter decreases for both peripheral and cen- 
tral collisions with a maximum at about 50% of the geo- 
metric cross section. Beyond 70%, the correction factor 
due to the reaction plane resolution is large, as shown the 
left-most panel in the Fig. 1. This limits the centrality 
range used in this analysis. 

The right-most panel in Fig. 1 shows the transverse 
momentum dependence of v-i for charged particles with 
respect to the reaction plane for minimum-bias events. 
The data above a px of 2 GeV/c clearly show a devi- 
ation from the monotonically increasing behavior seen 
at smaller px- The systematic errors are shown as line 
bands, which are estimated by several reaction plane 
methods using the two single beam counters or com- 
bined beam counters and by several different ways to 
correct non-uniform reaction plane distribution: "inverse 
weighting," "re-centering of sine and cosine summation," 
"Fourier expansion" and combinations of those above 
[11,18]. Those systematic errors are estimated to be 
about 10%, depending on centrality, and are indepen- 
dent of pt- Above 3 GeV/c, background tracks result in 
an additional systematic error of about 10%, depending 
on pt, which is included in the upper error band [19]. 

In Fig. 2, the transverse momentum dependence of 
V2 for identified particles is shown. The top-left panel 
shows negatively charged particles, while the top-right 
panel shows positively charged particles as described in 
the figure caption. The statistical errors and the system- 
atic errors are plotted independently. From the lambda 
particle spectra measured in the PHENIX central arm, 
it is determined that approximately 35% of the protons 
originate from lambda decays ( "lambda feed-down" ) [20] . 
The effect of the lambda feed-down on the measured vi 
of the proton is studied by varying the lambda i>2 with 
Monte Carlo simulation. Protons resulting from lambda 
feed-down increase the measured V2 value. Using the 
value of the lambda V2 measured at ^/s NN = 130 GeV at 
RHIC [3] , the effect on the proton V2 would be less than 
10%. Less than 5% of protons originate from decays of 
particles not involving the lambda. Based on further sim- 
ulations of their decays to protons, we estimate that the 
total systematic error due to feed-down is at most 11% 
depending on px, which is included in the lower sytem- 
atic error band in Fig. 2. 

The combined positive and negative particles are 
shown in the bottom-left panel. The lines in that panel 
represent a hydrodynamical calculation [21] including a 
first-order phase transition with a freeze-out tempera- 
ture of 120 MeV. The data show that at lower px (< 2 
GeV/c), the lighter mass particles have a larger v 2 at 
a given px, which is reproduced by the model calcula- 
tions. We note, however, that the difference between the 
charged kaons and charged pions is larger than the model 



predicts. 

A striking feature observed at higher px is that the 
v 2 of p and p are larger than for ir and K at p T > 
2 GeV/c. This is in sharp contrast to the hydrodynami- 
cal picture, which would predict the same mass-ordering 
for v 2 at all px- In our data the mesons begin to show 
a departure from the hydrodynamical prediction at px 
of about 1.5 GeV/c, while the (anti)baryons agree with 
the prediction up until 3 GeV/c but may be deviating 
at higher px. Such behavior is predicted by the quark 
coalescence mechanism [22], as shown in the bottom- 
right panel where both t> 2 and px have been scaled by 
the number of quarks. This could be an indication that 
the t>2 of measured hadrons is already established in a 
quark-matter phase, although it does not explain why 
the quark V2 would saturate with px- There exist other 
scenarios that could be applicable at RHIC, but we have 
selected two simple models (hydrodynamical and quark 
coalesence) only to emphasize the experimental evidence 
of the crossing of V2 for mesons and baryons. 

As an additional illustration of the different behavior 
for mesons and baryons, the transverse momentum de- 
pendences of the V2 parameter are shown in Fig. 3 for 
different particles and different centralities. Since the 
particle identification separation of K and p goes up to 
4 GeV/c, the combined 7r and K can be compared with 
protons up to 4 GeV/c. The charged particle acceptance 
is larger than the TOF acceptance where the particle 
identification can be performed. Therefore, the statis- 
tical fluctuations for the charged particle w 2 are smaller 
than for the p, p and tt + K. The trend exhibited in Fig. 2 
for minimum bias spectra, in which the i>2 parameter for 
(anti)baryons exceed those for mesons at px > 2 GeV/c, 
is shown here to occur for all centralities. 

In summary, the value of the V2 parameter for identi- 
fied and inclusive charged particle production at mid- 
rapidity has been measured with respect to the reac- 
tion plane defined in the forward and backward rapid- 
ity regions in ^/s NN = 200 GeV Au+Au collisions, us- 
ing the PHENIX experiment at RHIC. The value of v 2 
for charged particles decreases for both peripheral and 
central collisions with a maximum at about the 50 th 
percentile of the geometric cross section. We have ob- 
served that for charged particles V2 increases with px up 
to about 2 GeV/c, then starts to saturate or decrease 
slightly. However, the detailed behavior is different for 
different particle species. The lighter particles have larger 
V2 than the heavier particles for px below 2 GeV/c. This 
trend is partly reversed above 2 GeV/c where the proton 
and anti-proton have larger w 2 than mesons, a pattern 
which persists over all centralities. A hydrodynamical 
calculation can reproduce the mass-ordering and mag- 
nitude of V2 for the different particles in the region up 
to 2 GeV/c, but fails to reproduce either in the px re- 
gion above 2 GeV/c. As an alternative, we investigated 
the quark-coalescence scenario, in which the anisotropy 
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of the final-state hadrons is largely inherited from the 
anisotropy of quarks in a preceding quark-matter phase. 
The quark-coalescence model makes a definite prediction 
for a simple scaling behavior between the vi for mesons 
and for (anti)baryons, and this scaling behavior is largely, 
though not perfectly, borne out in our data. Further mea- 
surements extending to higher involving more identi- 
fied species will be required to discriminate among alter- 
native scenarios for the origin of elliptic flow at RH1C. 
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FIG. 1. (color online) Correlation of reaction planes be- 
tween two beam counters for the second moment is shown as 
a function of centrality (bottom-left) and the correction factor 
for the combined reaction plane resolution of two beam coun- 
ters is shown as a function of centrality (top-left). The value 
of V2 for charged particles is shown as a function of centrality 
(middle) and as a function of pt (right). 
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FIG. 2. (color online) Transverse momentum dependence 
of V2 for identified particles, ir~ , K~ , p (top-left) and tt + , 
K + , p (top-right). The circles show p and p, the squares show 
K + and K~ , and the triangles show ix + and 7r~ for minimum 
bias events. Statistical errors are represented by error bars 
and overall systematic error due to all sources by the solid 
lines in the top two panels. The combined positive particles 
and negative particles are shown in the bottom-left panel, 
and the lines there represent the result of a hydrodynamical 
calculation [21] including a first-order phase transition with a 
freeze out temperature of 120 MeV for n, K and p from upper 
to lower curves, respectively. The bottom-right panel shows 
the quark V2 as a function of the quark pt by scaling both 
axes with the number of quarks for each particle, as motivated 
by a quark coalescence model [22]. 
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FIG. 3. (color online) Transverse momentum dependence 
of V2 for combined ir~ and K~ (top) or 7r + and K + (bottom) 
compared with p (top) and p (bottom). In addition, results for 
inclusive negative (top) and positive (bottom) charged parti- 
cle distributions are plotted as open squares. From the left to 
right, the different centrality selections are shown for 0-20% 
(left), 20-40% (middle) and 40-60% (right). 
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